Qubit thermometry for micromechanical resonators 
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We address estimation of temperature for a micromechanical oscillator lying arbitrarily close to its quantum 
ground state. Motivated by recent experiments, we assume that the oscillator is coupled to a probe qubit via 
Jaynes-Cummings interaction and that the estimation of its effective temperature is achieved via quantum limited 
measurements on the qubit. We first consider the ideal unitary evolution in a noiseless environment and then 
take into account the noise due to non dissipative decoherence. We exploit local quantum estimation theory to 
assess and optimize the precision of estimation procedures based on the measurement of qubit population, and to 
compare their performances with the ultimate limit posed by quantum mechanics. In particular, we evaluate the 
Fisher information (FI) for population measurement, maximize its value over the possible qubit preparations and 
interaction times, and compare its behavior with that of the quantum Fisher information (QFI). We found that 
the FI for population measurement is equal to the QFI, i.e., population measurement is optimal, for a suitable 
initial preparation of the qubit and a predictable interaction time. The same configuration also corresponds to 
the maximum of the QFI itself. Our results indicate that the achievement of the ultimate bound to precision 
allowed by quantum mechanics is in the capabilities of the current technology. 

PACS numbers: 42.50. -p, 03.65.-w 



I. INTRODUCTION 

The edge between classical and quantum description of a 
phenomenon is related to the interactions occurring between 
the system under investigation and its environment. As a con- 
sequence, if we could, in ideal conditions, avoid irreversible 
interactions among them we should observe the emergence of 
quantum behavior even in macroscopic systems. As a mat- 
ter of fact, the technological developments of the recent years 
have made it possible to start inquiring into the quantum limit 
even in mesoscopic mechanical systems and experiments have 
been designed which realize a solid state analogue of cavity 
quantum electrodynamics. Many of these experiments focus 
on detecting the quantization of vibrational modes in a me- 
chanical oscillator fH- tlUl . Experimental conditions such that 
a mechanical object may behave in a quantum fashion are 
achieved in the low temperature regime. For example, for a 
single vibrational mode of energy haj to show quantum fea- 
tures, as the quantization of lattice vibrations, temperatures 
T <C jr^ are required, which for a micro-sized object oscillat- 
ing in the microwave band correspond to few mK. 

In this framework it has become increasingly relevant to 
have a precise determination of the temperature. However, for 
a quantum system in equilibrium with a thermal bath, there 
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is no linear operator that acts as an observable for tempera- 
ture. Temperature, thought as a macroscopic manifestation 
of random energy exchanges between particles, still retains 
its meaning but we have lost any operational definition. This 
kind of impediment often occurs in physics, and especially in 
quantum mechanics, whenever one is interested in quantities 
which are not directly accessible, i.e. they do not correspond 
to observable quantities. This may either be due to experimen- 
tal impossibilities, or be a matter of principle, as it happens 
for nonlinear functions of the density operator. In both cases, 
it turns out that the only way to gain some knowledge about 
the quantity of interest is to measure one or more proper ob- 
servables somehow related to the parameter we are interested 
in, and upon suitably processing the outcomes, to come back 
and infer its value. Hence, any conceivable strategy aimed to 
evaluate the quantity of interest ultimately reduces to a param- 
eter estimation problem. Relevant examples of this situation 
are given by estimation of the quantum phase of a harmonic 
oscillator 012M15I1 . the amount of entanglement of a bipartite 
quantum state itlq - UoTl an d the coupling constants of different 
kinds of interaction s lll9l429ll . Here we focus on the estima- 
tion of temperature 03011 and, motivated by recent experimen- 
tal achievements [ 1 1], we specifically refer to schemes where 
a micromechanical resonator is coupled to a superconducting 
qubit, and then a measurement of the excited state popula- 
tion is performed on the qubit itself. From the statistics of the 
population measurement is then possible to obtain informa- 
tion about the oscillator state, e.g. infer how close it is to the 
ground state, and in turn its temperature. 

In this context an optimization problem naturally arises, 
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aimed at finding the most efficient inference procedure lead- 
ing to minimum fluctuations in the temperature estimate. In 
this paper we address this problem in the framework of local 
quantum estimation theory (QET) IT32l - l37ll . We solve the dy- 
namics of the qubit-resonator coupled system and, in order to 
match realistic scenarios, we also take into account an effec- 
tive model for non dissipative decoherence. Then, we evaluate 
the Fisher information (FI) for the estimation of temperature 
via population measurement (hereafter referred to as the FI of 
the population measurement) and find both the optimal ini- 
tial qubit preparation and the smallest temperature value that 
can be discriminated. Moreover, we evaluate the Quantum 
Fisher Information (QFI) in terms of the symmetric logarith- 
mic derivative in order to calculate the ultimate bound to pre- 
cision allowed by quantum mechanics. This enable us to show 
that population measurement is indeed optimal for a suitable 
choice of the initial preparation of the qubit, and to provide 
quantum benchmarks for temperature estimation. 

It is worth noting at this point that we are not discussing 
here temperature fluctuations in a thermodynamical setting. 
Although temperature itself may not fluctuate, as it is sug- 
gested by quantum thermodynamical approaches 113 811 . we ex- 
pect that fluctuations always appear in the temperature esti- 
mates coming from indirect measurements [39, 40]. Quantum 
estimation theory provides the tools to evaluate lower bounds 
to the amount of fluctuations for a given measurement, as well 
as the ultimate bounds imposed by quantum mechanics. 

The paper is structured as follows. In Sec. [EI] we de- 
scribe the interaction model: first we briefly review the uni- 
tary Jaynes-Cummings dynamics for the coupled system and 
describe the measurements performed on the qubit, and then 
we take into account the decoherence effects. In Sec. Hill we 
show how QET techniques applies to our system, providing 
explicit formulas for both the FI and the QFI. The results are 
finally shown in detail in Sec. [IV]both for the unitary and the 
noisy dynamics. Sec. [V]closes the paper with some conclud- 
ing remarks. 



II. THE PHYSICAL MODEL 



field Hamiltonian reads: 



H F = h£l a) a 



(D 



where ft denotes the frequency of the vibrational mode. We 
assume the resonator in a thermal equilibrium state, i.e. de- 
scribed by the density operator 



Qf = 



exp(-/3ff F ) 
Tr [exp(-/3H,)] 



n=0 



p n {Sl,P)\n) (n\ 



where (3 = (fc^T)^ 1 and: 

p n (Q,0)=e-^ n 



1 



(2) 



The resonator is coupled to a superconducting qubit whose 
initial preparation is under control and, after a given interac- 
tion time, the excited state population is detected. The qubit is 
treated as a normalized vector in a two-dimensional complex 
Hilbert space Hq, with {|e) , \g)} providing an orthonormal 
basis. The qubit is initially prepared in a pure state 



e) + e^sin|| 5 ) 



(3) 



with ip € [0, 2ir) and £ [0, tt]. Hence the qubit density oper- 
ator reduces to the projector g Q = \ip) (ip\. Being a two-level 
system, by appropriately choosing the zero energy level and 
denoting by cj its transition frequency, the qubit Hamiltonian 
can be written as 

fiuj 

The qubit-resonator interaction is the interaction between a 
single-mode bosonic field and a two-level system. In the 
rotating-wave approximations and for the near-resonant case, 
i.e., for small values of the detuning S = us — SI we have the 
Jaynes-Cummings (JC) model with Hamiltonian 

H]c = H q + H F + H int 



a z + Mia} a + h\(a + a + er_at) . (4) 



As the temperature decreases a mechanical oscillator starts 
to exhibit its quantum nature, which mainly manifests itself 
in quantization of the vibrational modes. Hence, for our pur- 
poses the resonator can be regarded as a collection of phonons 
in a thermal equilibrium state. We assume that the resonator is 
built as to display an isolated mechanical mode at a given fre- 
quency, so that it can be modeled, rather than a phonon bath 
with some spectral distribution, as a single mode phonon field 
in thermal equilibrium. 



A. Unitary dynamics 



Let Hr be the infinite dimensional Hilbert space associ- 
ated with the single mode phonon field. Upon introducing the 
creation and annihilation operators [a, at] = 1 one has the 
number operator N — a^a, and its eigenstates {|n)}^_ . '^ rie 



The unperturbed Hamiltonian Hjq = H q + H F satisfies the 
eigenvalues equations 



H$ \k,n) = H 



\k,n) , 



with k — e.g and with the correspondences o e, 1 f-> g. 
In Eq. @)A £ M represents the coupling strength, a+a and 
cr-at stand respectively for the operators <r+ ® a, ct_ <g) 
acting on the tensor product space, where a± are the qubit 
ladder operators. Upon choosing a suitable rotating frame one 
rewrites the Hamiltonian in interaction picture Hjc- 



Hie 



+ h\ (cr+a + er_at) 



(5) 



The interaction only couples, for a given n, the states \e,n) 
and \g,n+ 1), and thus it is possible to study the interaction 
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inside the two-dimensional manifold spanned by these states 
leading to a representation - the so called dressed states basis 
- where Hjc is diagonal. We further assume the absence of 
any initial correlations between the qubit and the oscillator, 
thus choosing at time t = the following factorized density 
operator 

g(0) = g Q ® q f , 

whose dynamical evolution with respect to the JC Hamilto- 
nian is given by: 

Q(t) = U(t)Q(p)U*(t), 

withC7(t) =exp(-|.ffjct). 

Time evolution entangles the qubit and the resonator ll4lll 
and the probabilities for the qubit to be found in the ground or 
excited state are obtained via the Born rule as 



p(j\P) = Tr QF [g(t) \j){j\9h 



J = e, g 



(6) 



where p(j\f3) denotes the conditional probability of obtaining 
the value j when the value of the temperature parameter is /3. 
Upon introducing the following quantum operation: 



Qq 



g P = Tr F [U(t) g Q ®g F U\t)] 



(7) 



where £ : C(TLq) — > C("H Q ), Eq. (0 can be equally rewritten 
at the level of the qubit subsystem alone, namely: 



V(j\P) = Tr Q [Qp \j) 01] 



(8) 



In the following we will refer to g P as the probe state: It de- 
scribes the qubit subsystem at time t, obtained as the partial 
trace over the phonon field of the overall evolved state of the 
coupled system. Since it is a density operator on Hq it can be 
arranged in a 2 x2 density matrix. We have 



n=0 



Pn({l,P)(^ Be 3 
\ Uge Vgg 



where: 



cos 2 9 n t 
X 2 n 



A 52 ■ 2 



sin 2 6 n -it, 



1 



Qeg 



e~ lv sintf 



+ i- 



26 



1-1 1 



Qge 



cos 6 n t 
and Q gg 



.2* . „ 
1— sin ff n t 



1 - Qe 



(9a) 



(9b) 
(9c) 



with: 



B. Effects of decoherence 

A purely Hamiltonian dynamics doesn't match realistic fea- 
tures. In real-life scenarios quantum coherence is hard to 
achieve in mechanical objects, and can be maintained for rel- 
atively small times (« 10 _9 s ). Complete Rabi oscillations 
between the phonon and the qubit excitation involve only the 
first Rabi half periods, then a damping of the probabilities 
p(j\P) to I is observed: the most striking signature of de- 
coherence. Hence we include in our model the treatment of 
non dissipative decoherence occurring between the qubit and 
the resonator. Following Ref. [31] we consider an effective 
model provided by adding a power-law term in the thermal 
distribution, which leads to probe state matrix elements given 
by: 



Qij 



e T "*ftj 



1 



1 - e -T-*l 



being the matrix elements of Eq. (O, as evaluated for the 
unitary case, i, j G {e, g} and 



7ll = 6(1 + n) a . 



More explicitly 



1 X 

n=0 

9ge = Q* eg and g gg = 1 - g ee 



(Qe 



Qggi 



(10a) 

(10b) 
(10c) 



One can see that the dynamical evolution now drives the qubit 
towards the maximally mixed state, described by the density 
operator |. 



in. QUANTUM THERMOMETRY 

In this section we apply the tools of (local) quantum esti- 
mation theory (QET) to the coupled qubit-oscillator system. 
An estimation problem always consists in two steps: at first 
one has to choose a measurement and then, after collecting 
a sample of outcomes, one should find an estimator, i.e. a 
function to process data and to infer the value of the quan- 
tity of interest. In our case, temperature, expressed as /J, is 
the unknown parameter which has to be estimated from the 
sample of outcomes coming from measurements performed 
on the qubit. The results, a string of zeroes and ones for the 
case of population measurement, are distributed according to 
the probabilities p(j \0) = Qjj ofEqs. <[8J and© [orEq. (TTOt 
in presence of decoherence]. The Cramer-Rao inequality es- 
tablishes that the variance Var(/3) of any unbiased estimator is 
lower bounded by 



9 n (S,X) = -v^ 2 + 4A 2 (n + 1) 



Var(/3) > 



1 



MF{P) 



(11) 
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where M is the cardinality of the sample, i.e., the number 
of measurements, and F(j3) the so-called Fisher information 
(FI): 



[dpp(e\P)f [dp P (g\P)Y 
p(e\(3) p(g\f3) 



(12) 



Efficient estimators are those saturating the Cramer-Rao in- 
equality and their existence depends on the statistical model. 
However, independently of the statistical model we have that 
for sufficiently large samples, i.e., in the asymptotic regime 
M 3> 1, maximum likelihood estimators are always efficient. 

Quantum mechanically, the probability of obtaining the 
outcome j G {e,.g} from a measurement is given according 
tofheBornrulebypfj'l/?) = Tr [gpIX,], where the probe state 
q p = Q P (f3) parametrized by the unknown quantity /3 is re- 
ferred to as the quantum statistical model, and the collection 
of operators {IF,}, IF, > 0, J2j H = I is the probability 
operator- valued measure describing the measurement taking 
place on the qubit. In our case the qubit excited state popu- 
lation is probed and the measurement reduces to a projective 
one, |e) (e| and \g) (g\ = I — \e) (e\, i.e., we are measuring the 
Pauli operator u z = \e) (e\ - \g) (g\. 

Once the observable is fixed, we optimize the estimation 
procedure by maximizing the FI over the qubit state parame- 
ters, # and ip, as well as over the parameters driving the inter- 
action - i.e., the detuning S and the interaction time t. In other 
words, by employing the optimal qubit preparation and tun- 
ing the interaction parameters one may find a working regime 
achieving the maximum precision for that kind of measure- 
ment. 

On the other hand, one may also maximize the FI over all 
possible quantum measurements. Upon defining the symmet- 
ric logarithmic derivative (SLD) Lg as the selfadjoint operator 
satisfying the equation 



Lj3Q f 



= dag 



it is possible to show that the Fisher information F(/3) of 
any quantum measurement is upper bounded by the follow- 
ing quantity: 



F(J3) < G{fl) = Tr [gpLj] 



(14) 



which is called quantum Fisher information (QFI). QFI does 
not depend on the measurement carried on the qubit — indeed 
being obtained by maximizing over the possible measure- 
ment. It is rather an attribute of the family of states Q P (I3) 
parametrized by the temperature. Looking back to the 
Cramer-Rao inequality Eq.([m> one sees that QFI allows one 
to write its natural quantum version 



Var(/3) > 



1 



MGifi) 



(15) 



The above equation represents the Quantum Cramer-Rao 
bound (QCR), i.e. the ultimate bound to the precision al- 
lowed by quantum mechanics for a given statistical model 



Qp{P). An optimal measurement, i.e. a measurement whose 
FI F(/3) — G(/3) equals the QFI for the parameter /3, is given 
by the observable corresponding to the spectral measure of the 
SLD Lg. On the other hand, other kind of measurements may 
achieve optimality for the whole range of values of j3 or for 
a subset of values. Indeed, we will see in the following that 
population measurement is optimal for a suitable choice of the 
initial qubit preparation. We remind that for the estimation of 
a single parameter, as it is in our case, the QCR may be always 
attained, and an estimator saturating Ineq. ( [TBI is called effi- 
cient. The existence of an efficient estimator depends on the 
statistical model. However, independently of the statistical 
model, for sufficiently large samples, i.e., in the asymptotic 
regime M ^> 1, maximum likelihood and Bayesian estima- 
tors are always efficient. 

Upon diagonalizing the probe state one achieves the decom- 
position g P = g + (ip+\ + Q- IV*-) (VM an d i s a bl e to 
solve the equation for SLD 

{ip+\d B g P |V> + ) 
L P = \V+) {ip+\ 



Q+ 



Q- 



l^-XV-l 



[(ip + \dpg P \ip-)\ip+) (tp- 



g+ + q- 
+ (^_|c^p|v>+>hM<V+|], 

finally obtaining an explicit formula for the QFI 



(16) 



G{fi) 



Q- 



2k KV-|0/»V+>l +\(ip+\dpip-}\ 



(17) 



where 



\dp^±)=a p (e\^ ± )\e)+dp{g\i) ± )\g), 



(13) and 



{q+ 



g- 



(1 - 2g, 



Eq. dTTb contains a first term which resembles the FI and a 
second one, truly quantum in nature, which leads to the QCR 
and vanishes whenever \ip±) does not depend on f). 



IV. DYNAMICS OF THE FISHER INFORMATION AND 
OPTIMAL WORKING REGIMES 

In this section we report results for the qubit-resonator 
coupled system with physical parameters chosen in a range 
matching the experimental setup of Ref. Ill ill . More specifi- 
cally, we present a systematic study of the FI for population 
measurement as a function of the state and interaction parame- 
ters, carrying out numerical maximization and finding the op- 
timal working regimes. We also evaluate the QFI of the family 
of states Qp{0) and find the ultimate bound to precision, i.e. a 
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benchmark in order to assess the performances of qubit ther- 
mometry via population measurement. 

Hereafter we work with dimensionless quantities by rescal- 
ing times and frequencies in units of the coupling A. We thus 
substitute time, detuning and decoherence parameters by their 
rescaled counterparts 



t 



T = Xt, S i — >j = S/X, b i — >b = b/\. 



Effective detuning 7 will range in |-y | 6 [0, 1.5]. Also a di- 
mensionless effective temperature j3 is defined, provided by 
the substitution 

p 1 — > p ee pm . 

For convenience, we continue to term P and b respectively P 
and b. 



A. Resonant Hamiltonian regime 

Upon using the expression of the diagonal matrix elements 
in Eqs. (O we have evaluated the FI of Eq. dT2b . We start the 
discussion by considering the resonant case, i.e zero detun- 
ing, and analyze the effect of detuning afterward in this Sec- 
tion. For convenience we adopt the notation F(P) for the FI, 
but keep in mind the complete dependence F(P; r, 7) on 
both the qubit degrees of freedom and the parameters 7 and 
t which drive the coupling. Notice that F(P) does not de- 
pend on the qubit phase ip: its building-blocks are in fact the 
probabilities p(e\P) and p(g\P), whereas cp only appears in 
off-diagonal matrix elements. Varying the parameter § from 
7r to we span the entire class of qubit preparation, starting 
from |1), going trough a superposition and ending in |0). 

Let us now consider the system at a fixed value of the tem- 
perature, e.g. where the resonator is supposed to be very 
close to the ground state, say P = 10. The probabilities 
p(j\P) = Qjj evolve periodically in time according to Eq. ©, 
as the coupled system undergoes Rabi oscillations. The cor- 
responding behavior of the FI is shown in the upper panel of 
Fig. 03 The FI displays a robust maximum at the optimal time 
r max = ^ for ■d = 7T, corresponding to prepare the qubit in its 
ground state. This maximum is, at the same time, the global 
and the smoothest one. In fact, as soon as •& is moved from 
7r the FI suddenly drops to zero, except for a sharp peak cen- 
tered in r max , monotonically decreasing with respect to 1), as 
shown in the lower panel of Fig. Q] Another maximum of the 
same order of the global one can be found at r) = (0, 7f) 
but it is extremely peaked, thus representing a bad (unstable) 
choice for a possible measurement. Upon inspecting the tem- 
poral evolution of the excited state probability we found that 
p(e\P) has a minimum at r = r max , a fact which gives us 
a physical insight on the FI behavior: since our goal is the 
estimation of a vanishing quantity which carries information 
about thermal disorder, we expect to find the maximum sen- 
sitivity in our predictions where the excitation is most likely 
stored - as a phonon - in the resonator, i.e., when p(e\p) is 
minimum. 
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FIG. 1: (Color online) Upper panel: FI for (3 = 10 as a function 
of the effective time r, for different $ values: d — tt (dashed blue), 
$ = 0.95 7r (dot-dashed magenta) and j? = (solid green). FI takes a 
pronounced global maximum at (#, r) = (jr, ?) while it is possible 
to see a secondary extremely peaked maximum, which occurs for 
t = 7r and preparing the qubit in |0). Lower panel: log-linear plot 
of the FI for j3 = 10 as a function of 1? for, r = ^ (dashed blue), 
T = 5 + e (dot-dashed magenta), r = it (solid green), r = tt + e 
(dotted red), with e = 0.01. 



1 

0.001 

10~ 6 
10~ 9 

10- 12 




FIG. 2: (Color online) Log-linear plot of the FI as a function of ef- 
fective time r for different values of /3. The qubit is prepared in the 
ground state |1) (1? = tt). From bottom to top ft = 15 (solid blue), 
P — 10 (dashed magenta), p = 5 (dot-dashed green), P — 1 (dotted 
red). Upon raising the temperature the FI no longer keeps a scale-free 
shape: thermal excitations modifies its profile making it irregular. In 
particular the global maximum comes earlier in time. 
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Let us now turn our attention to the dependence of the FI 
on the temperature itself. In Fig.|2]we show, on a logarithmic 
scale, the temporal evolution of the FI for different values of 
8. FI varies over several orders of magnitude, matching our 
intuition that the closer we are to the ground state, the harder is 
to achieve a given precision in estimation of temperature. Fur- 
thermore, upon lowering the temperature, the temporal evolu- 
tion of p(j\f3) becomes less involved, finally approaching the 
exactly periodic one of Rabi oscillations, which in turn freezes 
the profile of the FI in a shape independent on the temperature 
itself. 

The qubit preparation 8 = ir is universally optimal, i.e., it 
leads to a maximum of the FI independently of the interaction 
time. After fixing = ir we have numerically maximized 
F(j3i) with respect to r. The solid blue line of the upper panel 
of Fig.|3]is the the log-plot of 

F Af ( / 9)=maxF(/9), 

T 

as a function of 8, from which it is apparent the exponen- 
tial decrease of the maximum value achieved by the FI for 
increasing j3. The Cramer-Rao inequality immediately relates 
this fact to an exponential loss of sensitivity moving towards 
the quantum ground state of the resonator. An other interest- 
ing feature that emerges from the maximization is a shift in 
the value of the optimal interaction time. In the lower panel of 
Fig. [3] we can recognize the existence of a steady value for the 
optimal time T max = 5 when approaching the ground state, 
while for smaller values of j8 the optimal time comes earlier. 
In fact, the temporal evolution of FI (see Fig. I3 not only pre- 
dicts an exponential increase of the global maximum when 
temperatures are raised, but also a shift of its location. 
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FIG. 3: (Color online) Upper panel: log-log plot of the FI maximized 
over t as a function of /?, with 9 — tt for different values of detuning: 
7 = (solid blue), 7 = 1 (dashed magenta), 7 = 1.5 (dot-dashed 
green). Bottom panel: the times r max which maximizes the FI as a 
function of /3, with 8 = n for different values of 7 (same values and 
colors of the upper panel). 



B. Effects of detuning 

In this section we take into account the possible existence 
of a nonzero detuning 7 between the oscillator and the qubit 
frequencies. This has two main consequences, which are both 
illustrated in Fig. [3] On the one hand, the maximum achiev- 
able value of the FI slightly decreases and, on the other hand, 
the optimal interaction time r max at which the maximum takes 
place anticipates. Therefore, the best working conditions to 
achieve the optimal sensitivity in the estimation of 8 corre- 
spond to have the qubit and the resonator in resonance. It is 
also worth to notice that 7 does not represent a critical pa- 
rameter, as the initial preparation of the qubit, since the FI 
dependence on 7 is smooth. One can see this in the upper 
panel of Fig. [3] where we see that curves corresponding to 
quite different values of the detuning are almost superposed. 

C. Quantum Fisher information 

In order to assess the performances of the population mea- 
surement in the estimation of temperature we have evaluated 
the QFI of the family g P {0). The diagonalization of the probe 
state has to be carried out numerically, hence in general an- 
alytical expressions of the QFI are not available. A first fact 



is that G(B) turns out to be independent on the qubit phase 
ip, which then does not represent an extra degree of freedom 
whereby gain more restrictive bounds to precision on Var(/3). 
Even the optimal qubit preparation for to the best conceivable 
measurement involves control of the parameter only. 

As we have done for the FT, we start to inspect the QFI be- 
havior for a fixed value of temperature 8 in the resonant case. 
Also for the QFI the maximum is achieved by preparing the 
qubit in the state \g) and probing it at time r max . In this case 
the behavior of G(8) is identical to that of F{8), as it is ap- 
parent by comparing Figs. Q~|and[4] In other words, for a given 
value of the parameter j8 into the range explored, the choice 
t) = (tt, r max ) makes population measurement optimal. 
Moreover, the QFI itself reaches its global maximum for that 
choice. Thus, provided that an optimal estimator is employed, 
e.g. maximum likelihood in the asymptotic regime, this strat- 
egy provides optimality in sense that either inequality (TT~4b is 
saturated and the right-hand side of QCR is as low as possible. 

This conclusion is confirmed upon a closer inspection of 
the probe state. When $ = tt the off-diagonal terms vanish 
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and g P is diagonal, with eigenvalues 



Q+ = Vn (0, P) sin 2 V7 2 + 4n 

n=0 
= 1 - Q + 

As a consequence, the QFI reduces to 



n + 7 2 /4 



(18a) 
(18b) 



£4 



Q- 



which coincides with the FI ruling the estimation of f3 via pop- 
ulation measurement. 

On the other hand, some striking difference emerges be- 
tween the performances of population measurement and that 
of the optimal one if the qubit is not prepared in the optimal 
(ground) state. 
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FIG. 4: (Color online) Upper panel : QFI for /3 = 10 as a function 
of r, for 1? = tt (dashed blue), $ = 0.95 tt (dot-dashed magenta) 
and $ = (solid green). QFI behaves like FI for -3 — tt leading 
to the same maximum, while for smaller angles it shows a smoother 
profile. For angles < $ < tt one may find measurements which 
improve the precision of temperature estimation. Bottom panel: QFI 
for f3 = 10 as a function of -d for r = ? (dashed blue), r = ? + e 



(dot-dashed magenta), r 
with e = 0.01. 



tt (solid green), r = tt + e (dotted red), 



In the lower panel of Fig.|4]we show G(f3) as a function of 
t for different values of fr. for § < ir the decrease of G is def- 
initely smoother than that of F and thus, in principle, some 
measurement may be found making the initial preparation a 



less critical parameter. Moreover inspecting the cut of the QFI 
along r = tt we note that the maximum in $ = becomes 
more achievable compared to the one of F(/3). All these fea- 
tures suggest that for qubit preparations different from the 
ground state there will be a sensible difference between the 
precision provided by population measurement and the opti- 
mal one implementable on the system. On the other hand, 
being the overall maximum achievable with population mea- 
surement, our results indicate that the achievement of the ulti- 
mate bound to precision allowed by quantum mechanics is in 
the capabilities of the current technology. 



D. Effects of decoherence 

In this section we discuss the solution of the reduced qubit 
dynamics in the presence of dissipative decoherence, see Eq. 
([Tol l, and inspect the corresponding behavior of the FI. For 
the sake of simplicity we consider zero detuning. Analogue 
results are obtained when including the detuning. 

The probabilities p(j\f3) — Qjj are damped so that, waiting 
for a sufficient long time, whose value depends on a and b, 
we would find them to be identically 1/2 or, equally stated, 
the dynamical evolution brings the state to the maximally 
mixed one. The contribution of decoherence is of the kind 
exp[— 6(1 + n) a r] for every n, where b has been rescaled in 
coupling units b 1 — > b/X. Being a multiplicative coefficient, 
as soon as b is different from zero, the exponential term will 
participate in killing the sums. Our calculations show a rel- 
evant dependence of the FI on the parameter b, namely val- 
ues b rs 10~ 5 are sufficient to produce visible effects, while 
varying a in the range (0,1) does not deeply influence of FI 
behavior. 

In Fig. [5] we show the temporal evolution of the FI for 
(3 = 10, in the presence of decoherence and for different ini- 
tial preparations of the qubit. In the Hamiltonian regime for 
large (3 the resonator is close to the ground state, the evolution 
of p(j\/3) is periodic and hence, due to Eq. ( fT2l . the same is 
true for the FI. Upon incorporating decoherence we see that 
FI decays at a rate depending on b and thus an irreversible 
dynamics emerges, which matches the physical evidence of a 
limited coherence time. On the other hand, a clear maximum 
at t = 7r/2 still appears, with a slightly decreased value of 
F(j3). In the lower panel of Fig. [5] we show the maximum 
value Fm = max T F{j3) for different values of the decoher- 
ence parameter. As it is apparent from the (log-log) plot for 
high temperature (smaller /3) the effect of decoherence is neg- 
ligible, whereas for increasing (3 the effect is becoming more 
and more relevant. 



V. CONCLUSIONS 

The temperature of a physical object cannot be directly 
measurable. On the other hand is can be regarded as a pa- 
rameter whose value can be indirectly inferred by measuring 
some proper observable and then suitably processing the out- 
comes, an inference procedure usually referred to as an esti- 



8 



10 
0.4 



0.3 
0.2 
0.1 







/ >» 




' * 1 




/ !! 




! 


i i 


■I 


i ii '. 




i ; 






i 


i « 




/ !! \ 


/ \ 


i hi 




i : 
i ! 


i i 


/ if. \ 


i i 


i 

i i 


! i \ 


i 






» i 


i \ 




1 ! 
i 

—' 


\ \ 


' > 1 1 



I I 

i i 
/ i 
/ i 



0.001 




FIG. 5: (Color online) Upper panel: Fisher information F(/3) for 
/3 = 10 as a function of r in the presence of decoherence and for 
different qubit preparations. The decoherence parameters are chosen 
as to a = 0.1 and b = 10 -5 . Dashed blue line stands for # = 7r, dot- 
dashed magenta for i? = 0.95 n while solid green ones for # = 0. 
Having included decoherence treatment enables us not to restrict the 
evolution to the first Rabi half-period. Lower panel: log-log plot of 
the Fisher information Fm (/3) maximized over the interaction time, 
and in the presence of decoherence, as a function of /3 and for fixed 
# = 7T, for 6 = (solid green), b — 10 -5 (dashed magenta), b = 
10~ 4 (dot-dashed blue). 



mation procedure. In the case of a micromechanical oscilla- 
tor with an isolated vibrational mode, effective schemes have 
been suggested and realized ll ill which rely on coupling the 
resonator to a superconducting qubit and probing the latter 
using population measurements. In other words, the qubit is 
employed as a quantum thermometer to demonstrate that the 
resonator has been cooled to its quantum ground state. In this 
paper we have analyzed in details qubit thermometry in these 
systems, i.e., the estimation of temperature via quantum lim- 
ited measurements performed on the qubit. In the framework 
of quantum estimation theory we have analyzed precision as 
a function of both the qubit initial preparation and the interac- 
tion parameters, and we have evaluated the limits to precision 
posed by quantum mechanics to qubit thermometry. 



We have computed the FI for population measurement, 
which is the appropriate figure of merit to assess the preci- 
sion of estimation, and have found that its maximum, and 
hence the minimum variance in the estimated temperature, is 
achieved by preparing the qubit in the ground state, and prob- 
ing it at an emergent time r max , which is predictable. Further- 
more, we have analyzed in details how the maximum depends 
on the temperature itself, on the detuning, and on the noise 
parameter when one takes into account non dissipative deco- 
herence. In order to evaluate the ultimate bound allowed by 
quantum mechanics to the sensitivity of temperature estima- 
tion, we have also computed the quantum Fisher information. 
We found that QFI is maximized for the same choice of qubit 
preparation and measurement time of the FI, and that for these 
common values the maxima of FI and QFI coincide. We thus 
conclude that population measurement is optimal for temper- 
ature estimation. 

The range of parameters addressed in our analysis is that 
of recent experimental implementations 111 ill . We thus con- 
clude that optimal estimation of temperature can be done with 
current technology. Since the FI of population measurement, 
and the QFI of the model, both decrease with the decrease 
of temperature, the estimation of lower temperature will be 
intrinsically less precise. On the other hand, since the are 
regimes, also in the presence of decoherence, where the max- 
ima of the FI and the QFI are reasonably smooth as a function 
of the qubit preparation and of the interaction time we do not 
expect any "no-go" theorem for temperature estimation. In 
other words, we expect that optimal estimation of lower res- 
onator temperatures, perhaps achievable with further experi- 
mental advances, will be still possible with population mea- 
surements. On the other hand, "optimality" will correspond 
to an inherently less precise procedure compared to the case 
of higher temperature. 

Our analysis shows the optimality of feasible qubit ther- 
mometry in providing quantum benchmarks for high precision 
temperature measurement, as well as an efficient operational 
quantification of temperature for mechanical modes lying ar- 
bitrary close to their ground state. In other words, achieve- 
ment of the ultimate bound to precision allowed by quantum 
mechanics is in the capabilities of the current technology. Our 
results also confirm that QET is a useful tool for assessing 
and comparing inference procedures arising in quantum lim- 
ited measurements M42I1 . even when mesoscopic objects are 
involved. 
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